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ABSTRACT 

We studied the clustering properties and multiwavelength spectral energy distributions of a complete sam- 
ple of 162 Lya-emitting (LAE) galaxies at z ~ 3.1 discovered in deep narrow -band MUSYC imaging of the 
Extended Chandra Deep Field South. LAEs were selected to have observed frame equivalent widths >80A 
and emission line fluxes > 1.5 x 10~'^ergs cm"^ s~'. Only 1% of our LAE sample appears to host AGN. The 
LAEs exhibit a moderate spatial correlation length of ro = 3.6;!;j' gMpc, corresponding to a bias factor b=l .T^'l, 
which implies median dark matter halo masses of logjoMmed = 10.9;!;o ^Mq. Comparing the number density of 
LAEs, 1.5 ± 0.3 X 10~^Mpc"^, with the number density of these halos finds a mean halo occupation ~1-10%. 
The evolution of galaxy bias with redshift implies that most z = 3.1 LAEs evolve into present-day galaxies with 
L < 2.5L*, whereas other z > 3 galaxy populations typically evolve into more massive galaxies. Halo merger 
trees show that z = Q descendants occupy halos with a wide range of masses, with a median descendant mass 
close to that of L*. Only 30% of LAEs have sufficient stellar mass (>^ 3 x lO^M©) to yield detections in 
deep Spitzer-IRAC imaging. A two-population SED fit to the stacked f/By/?/z7/rH-[3.6,4.5,5.6,8.0]^m fluxes 
of the IRAC-undetected objects finds that the typical LAE has low stellar mass (1.0^o 4 x lO'^Mo), moderate 
star formation rate (2 ± IMoyr"'), a young component age of 20;!;[QMyr, and little dust (Ay < 0.2). The best fit 
model has 20% of the mass in the young stellar component, but models without evolved stars are also allowed. 

Subject headings: galaxies:high-redshift - galaxies:formation - galaxies:evolution - large-scale structure of uni- 
verse 



1. INTRODUCTION 

The discovery of high-redshift Lya-emitt ing galaxies 
(LAE s ) opened a new frontier in astronomy dCowie & Hul 
119981: IHu et all 11998 !). Because the Lya fine is easily 
quenched, a galaxy with detectable Lya emission is likely 
dust-free, i.e., in the initial phases of a burst of star formation. 
The Lya lines have large equivalent widths (20A<EWiest <^ 
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lOOA) and broad velocity widths (100 km s"' <FWHM<800 
km s"') and are often asymmetric, indicative of high- 
redshift galaxies un dergoing active star-forma tion (e.g.. 
Manning et all 120 00^, Kudritz ki et all 12000; Arnaboldi et al.l 



20021: iRhoadset al. 2003; Daw son et al.ll2004l:|Hu et al.ll2004t 
Vene mans et al.1 l2005i: ,Matsuda et all 120061: bronwall et all 



2007). 

Other high-redshift galaxy populations (including Lyman 
break galaxies. Distant Red Galaxies, Sub-Millimeter Galax- 
ies) exhibit strong clustering and sh ould evolve into elliptical 
and giant el liptical galaxies today jAdelberger et al.l l2005bl: 
lOuadri etal . 2007; Webb et al. 2003]). These objects were se- 
lected by unusually strong rest-frame continuum emission in 
the ultraviolet, optical, an d far-IR respective ly, resulting in 
IO'^Lq < Lh„i < IO'^Lq jReddv etalj|2006h . Such sti-ong 
continua appear to occur primarily in deep potential wells that 
are strongly biased versus the general distribution of dark mat- 
ter halos. LAEs instead offer the chance to probe the faint end 
of the (bolometric) luminosity function of high-redshift galax- 
ies, which contains the majority of galaxies. The strong Lya 
emission line allows detection and spectroscopic confirmation 
of LAEs with typical bolometric luminosities ~ lO'^L©. A 
detailed calculat ion of the LAE lu minosity function at z ~ 3.1 
IS given in Gronwalletal.1 (|2007j). 

Spectral energy distribution (SED) modelling of the stacked 
UBVRIzJK photometry of 18 LAEs in the Extended Chandra 
Deep Field-South (ECDF-S) (iGawiser et alJl2006ah showed 
the average galaxy to have low stellar mass (< lO' Mq) and 
minimal dust (see also Nilsson et al. 2007). LAEs have the 
highest specific star formation rates (defined as SFR divided 
by st ellar mass) of an y type of galaxy, implying the youngest 
ages jCastro Ceron et al. 2006). Because Lya emission is eas- 
ily quenched by dust, LAEs have often been characterized 
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as proto galaxies experiencing th eir first burst of star forma- 
tion (e.g. lHu & McMahonll 19961) . However, the differing be- 
havior of Lya and continuum photons encountering dust and 
neutral gas makes it possible for older galaxies to exhibit 
strong Lya emission when morphology and kinematics fa- 
vor the escape of these photons tow ards Earth ( Neufeld 19911 
iHaiman & Spaansll 19991: [Hansen & Oh 2006). This could al- 
low older, dusty galaxies with actively star-forming regions to 
exhibit Lya emission with high equivalent width. 

SED modelling of LAEs using Spitzer-IRAC dFazio et alj 
12004) to probe rest-frame near-infrared wavelengths, where 
old stars dominate the emission, has yielded mixed results. 
iPirzkal et al.l (l2007h report extremely young ages of a few 
Myr and low stellar masses (lO'^M© <M* < lO^M©) from 
SED modelUng of 9 LAEs at 4.0 < z < 5.7 i n the Hubble Ul- 
tra Deep Field. However, iLai^J: alj i2007a| ) performed SED 
fitting to 3 LAEs with IRAC detections out of a sample of 12 
z = 5.7 LAEs in GOODS-N and found ages as high as 700 
Myr and significant stellar masses (10^ < M0 < 10'"), mak- 
ing it appear that these LAEs were not undergoing their first 
burst of star formation. The 9/12 LAEs lacking IRAC detec- 
tions are presumably younger and less massive and might lack 
an evolved population. Investigating the nature of the LAEs 
without IRAC detections requires a stacking approach to see 
if the typical LAE stellar mass is low enough to have been 
generated in a single ongoing starburst. Stacking will yield 
the best results when applied to a large statistical sample of 
L AEs in a region wi th dee p IRAC imaging. 

iGronwall et alJ (l2007h present the largest avail- 
able sample of LAEs in an unbiased field, 162 
LAEs at z = 3.1 in the ECDF-S discov e red as 
part of the MUSYC survey (Gawis er et all l2006bl 
[FEtpT77www.astro.yale.edu/MUSYC'). The ECDF- 
S has been targeted with deep narrow-band imaging and 
multi-object spectroscopy, complemented by public broad- 
band UBVRIzJK, Spitzer-HlRAC and ChandraH-AClS-1 imag- 
ing. We improve the constraints of iGawiser et al.l (l2006ah 
using the MUSYC UBVRIzJK photometry of our larger sam- 
ple of LAEs and adding IRAC [3.6^m,4.5/im,5.8/im,8.0^m] 
Cycle 2 legacy images from SIMPLE (Spitzer 
IRAC/MUSYC Public Legacy of the ECDF-S, 
|http : / /w ww . astro . yale . edu/ dokkum/ simple). 

This Letter summarizes our imaging and spectroscopic ob- 
servations of LAEs, presents our results from clustering anal- 
ysis and SED modelling, and discusses the implications for 
the formation process of typical present-day galaxies. We 
assume a ACDM co smology consistent with WMAP results 
(lSpergeletalJl2007h with n,„ = 0.3, f^A = 0.7, Hq = 70 km s"' 
Mpc , and rms dark matter fluctuations on 8/1"' Mpc scales 
given by erg = 0.9 . All correlation lengths and number den- 
sities are comoving. We have suppressed factors of /170 in 
reporting correlation lengths, number densities, dark matter 
masses, stellar masses and star formation rates. 

2. OBSERVATIONS 

Our narrow-band 4990A and U BVRIzJK im- 
ages of ECDF-S are described in 'Gronw all et alJ 
(l2007h and iGawiser et alj (|2006a) and are available at 
|http : / /www, astro, yal e . edu/MUSYC The final 
images cover 31.5' x 31.5' = 992 arcmin"^ to a narrow -band 
completeness limit of ^ 1.5 x 10"'^ ergs cm~^ s~' (AB=25.4 
in the 50A FWHM NB4990A filter). Figure [U shows our 
complete sample of 162 strong Lya emitting galaxies at 
z ~ 3.1 with equivalent width >80A in the observed frame. 



28 LAEs lie in the region surv eyed by the GOODS Legacy 
program (iDickinson et al.ll2003h . 
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Fig. 1 .— Plot of the 31.5' X 31.5' Extended Chandra Deep Field-South, 
showing the MUSYC BVR-selected catalog of 84490 objects as tiny dots. 
The rectangle shows the GOODS-S region, and the inner square shows the 
Hubble Ultra Deep Field. Lya-emitting galaxies at z = 3. 1 are shown as solid 
circles. 

Multi-object spectroscopy of 92 LAE candidates, along 
with other MUSYC targets, was performed with Magellan- 
Baade-HlMACS on Oct. 26-27, 2003, Oct. 7-8, 2004, Feb. 
4-7, 2005, Nov. 2-3, 2005, Oct. 25-27, 2006, Nov. 21-22, 
2006 and Feb. 18-20, 2007. The 300 line/mm grism was 
used with 1.2" sHtlets to cover 4000-9000A at a resolution 
of R = 640 i.e., 470 km s"', at the wavelength of Lya emis- 
sion. Mask exposure times ranged from 2 to 5 hours, with the 
longer exposures sufficient to detect Lya emission lines down 
to our completeness limit of 1.5 x 10"'^ ergs cm~^ s~', as- 
suming clear conditions and minimal slit losses. Details of 
our spectroscopy will be given in P. Lira et al. (in prep). Red- 
shifts were confirmed to lie at 3.08 < z < 3.14 for 61 of the 
LAEs, with 1 interloping AGN at z = 1 .60 where [C 111] A 1909 
falls in the narrowband filter, and the other 30 objects lacking 
sufficient S/N to yield redshifts. Our success rate for the slit- 
masks with the highest S/N was 90%, setting an upper limit 
of 10% on possible contamination of our LAE sample. The 
rate of non-detections was higher in masks with shorter expo- 
sure times resulting from weather or instrument challenges, 
consistent with the reduced S/N. Our spectroscopy shows that 
the sample is not contaminated by z = 0.34 [O 11] emission- 
line galaxies, which are the typical interlopers for narrow- 
band-selected LAE samples. These have been eliminated by 
requiring observer's-frame EW> 80A which ehminates all 
but the rarest [0 11] emitters jTerlevich et alJl99lMHogg et alj 
IT998llStern et alJi2000 ). 

The Lya emission in LAEs appears to derive from star for- 
mation rather than AGN activity; only 2/162 objects in our 
complete sample are detected as X-ra y sources in the Chan- 
dra catalogs of CDF-S and ECDF-S jAlexander et"an 120031: 
lLehmeretalll2005l: IVirani et alj[2006h . One is the z = 1 .6 in- 
terloper but the other is at z = 3.092. One additional object at 
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z = 3.1 17 is detected in X-ray photometry at the narrow -band 
source position. At z = 3.1, Chandra receives X-ray emission 
from 2-30 keV (rest-frame), meaning that Compton thick ob- 
scuration (A^H >'~-" lO^-'cm"^) is needed to hide AGN. Even 
such heavily obscured AGN are likely to reveal their presence 
via narrow emission lines, which should be indicated by high- 
ionization UV lines like C IV accompanying Ly a. Amongst 
our LAE spectra, only the three X-ray sources show signs of 
AGN activity in the form of emission lines other than Ly a, 
and th e other 159 LAEs are undetected in a stacked X-ray 
image jGronwall et al.ll2007 h. We therefore expect that very 
few LAEs contain AGN that dominate their Lya or contin- 
uum emission. Our two X-ray sources at z = 3.1 imply that 
AGN are present in only 1 .2 ± 0.9% of Ly a selected galaxies 
at this redshift. We restrict our subsequent analysis to the 159 
LAEs without X-ray detections. 

3. CLUSTERING ANALYSIS 

We used the iLandv & Szalavl ( Il993h estimator to measure 
the angular correlation function using histograms of pairs of 
points at separation 6 between the data catalog (D) and itself 
{DD) as well as the cross-correlation and autocorrelation with 
a set of random catalogs {R). We modelled this observed cor- 
relation function as an intrinsic 1^(6*) minus the "integral con- 
straint" caused by estimating the sky density of LAEs from 
our own survey (see .Peeblea, 1980; infante.. 1 994„ H. Francke 
et al. in prep) 



DD{9)-2DR(9) + RR{e) 
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where a"^ = — j j w(9i2)d^idfl2 
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Fig. 2. — Angular auto-correlation function data with best-fit model (solid 
line) and best-fit spatial congelation length and dark matter halo bias and mass 
parameters listed. The bin below 30" (square) is not used in the fit. 

Figure |2] shows our binned data along with the best- 
fit model. Error bars show the uncertainties in each 
bin of a„(9) = (1 + w(9)) / ^/RR(9) (iLandv & Szalavl [T99l 
iGawiser et al.ll20()6bh . The analysis was restricted to scales 
larger than 30" («1 comoving Mpc) that are insensitive to 
the possible presence of multiple LAEs in some dark matter 
halos. In order to determine the expected redshift distribu- 
tion, Nexpiz), we performed a Monte Carlo simulation plac- 
ing a large number of LAEs over the redshift range 3 < z < 



3.2 and assigning them equivalent widths drawn at random 
fror n the equivalent width distribution observed for our sam- 
ple (lGronwalletal.1120071) . We then used the NB4990 filter 
bandpass to calculate the excess narrow-band flux that would 
be observed, removed all LAEs with "observed" equivalent 
width <80A, and measured the redshift distribution of the se- 
lected objects. Fig. [3]shows that the expected redshift distri- 
bution is narrower than the transmission curve. This occurs 
because LAEs far from the central redshift must have very 
high equivalent widths to be selected through narrowband ex- 
cess, due to the reduced filter transmission at the wavelength 
of their Ly a emission lines. Fig. |3]shows that the histogram 
of observed LAE redshifts is consistent with Nexpiz)- The only 
bin inconsistent with poisson fluctuations of the expected red- 
shift distribution at the 9is at 3.085 < z < 3.090, which dis- 
plays a 3cr excess, revealing an overdensity at this redshift. 
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Fig. 3. — Histogram of LAE redshifts in ECDF-S. The thin solid curve is 
the filter transmission curve normalized to have the same ai'ea, and the thick 
solid curve is the expected redshift distribution described in the text. 



We used N^r piz) to deproje ct the angular correlation func- 
tion, following Simon] (|2007|) . In order to avoid degeneracy 
between the clustering length ro and the power-law index 7 
in the underlying spatial correlation function, ^(r) = (r/rQ)~^, 
we assumed a typical power-law with 7 = 1.8. This yielded a 
moderate clustering length, ro = 3.6;!;° f, Mpc (comoving). The 
narrow redshift distribution of narrow-band selected LAEs 
reduces the loss of angular clustering signal due to projec- 
tion, allowing a high S/N measurem ent of moderate cluster- 
ing. Following lOuadri et al.l (12007) . this value of ro corre- 
sponds to stronger clustering than the dark matter at z = 3 . 1 
by a bias factor b=l ■7'!^)l- Note that bias factors are robust to 
the degeneracy between ro and 7. This bias factor is shared by 
the population of dark matter halos with mass es greater than 
" Mq (ISheth & Tormen|[T999i) . implying a 



logioMmin 



1 



median dark matter halo mass of logjgMnied 



10.9t5 Mq. 

If 10% of the LAEs were unclustered low -redshift contami- 
nants, the corrected value of ro would be 10% higher, yielding 
halo masses ^ 50% higher. 
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The comoving n umber density of our L AE sample is 1 .5 ± 
0.3 X 10"^Mpc"3 (iGronwall et alJ l2007h . where the uncer- 
tainty includes variance due to large-s cale structure in our sur- 
vey volume for objects with b = l.l dSomerville e t al. 2004]). 
The number density of the corresponding dark matter halos 
is 3tl^ X lO'^Mpc"-*, implying a "mean halo occupation" of 
5^4"% for the LAEs. There is significant freedom in how 
the LAEs could be assigned to this subset of the available 
dark matter halos. However, the LAE median halo mass must 
roughly follow the result logjQ Mmed = lO.Q^g 9 Mq in order to 
explain the observed clustering bias.'^ Explorations of com- 
plex halo occupation distribution (HOD) models show that the 
assumption of one galaxy per halo made in our determination 
of Mined can ca use additional uncertainties of up t o 0.2dex at 
2 = 3.1,/?= 1.7 (iLee et al.ll2006tlZheng et alJl2007h . 

4. SED MODELLING 

iLai et alj (l2007bh offer a detailed description of our IRAC 
photometry along with single-component SED fitting of the 
detected and undetected objects and a comparison of their 
continuum properties with those of Lyman break galaxies. 
76 of our LAEs fall within regions of the SIMPLE images 
(which include the GOODS IRAC images) where the lack of 
bright neighbors enables IRAC photometry accurate to very 
low fluxes. Only 24 LAEs (30%) are detected by IRAC at 
fluxes above the 2a SIMPLE flux limit; these objects repre- 
sent the high-mass end of the LAE mass function and appear 
to have stellar masses > 3 x lO'^M© Only 2 of these LAEs are 
detected in our J,K images, which are two magnitudes shal- 
lower than the IRAC 3.6,4.5 /im images. The IRAC-detected 
LAEs are brighter in the rest-UV and rest-optical continuua, 
with mean R-band and 3.6/im fluxes corresponding to magni- 
tude 25.4 and 24.4 respectively, compared with 26.7 and 26.6 
for LAEs not detected by IRAC. In order to investigate the 
full SED of typical LAEs, which are too dim to be detected 
individually in our NIR and Spitzer images, we measured av- 
erage fluxes from stacked images of the 52 LAEs (70%) lack- 
ing IRAC detections. We show the resulting SED in Fig. 21 
where the V-band flux has been corrected for the contribution 
of the Lya emission lines to this filter Uncertainties in the 
stacked photometry were determined using bootstrap resam- 
pling to account for both sample variance and photometric 
errors. 

Instead of modelling LAEs with a single stellar population, 
we analyzed the extent to which the data allow the presence 
of an underlying evolved population. We adapted the method 
of K. Schawinski et al. (in prep.) to model the star forma- 
tion histories using a two-burst scenario with the old compo- 
nent as an instantaneous burst and the young component as 
an ex ponentially dec lining starburst with variable e-folding 
time. iMarastonI (12005 ) population synthesis models were used 
with m etalhcity ranging from 0.02 sol ar to solar. a | Salpeteij 
( Il955h initial mass function, and the Calzetti et alj ( l200d) 
dust law. The best-fit model shown in Fig. |4] corresponds 
to a stellar mass of 1.0;!;o4 ^ ^0' ^'■^ formation rate of 
2 ± lM0yr"\ and dust extinction Ay = O.OiliQ g (only positive 
values of Ay were considered). Figure |5] shows the results 
for the age of the young population versus the mass fraction 
of the young population. The young population has an age 
of 20^JoMyr with an e-folding timescale t = 750 ± 250Myr 

The quantity that must be preserved is the mean halo bias. The median 
halo mass is a simpler statistic that is also robust in typical HOD models, and 
the difference is far smaller than the reported uncertainties. 
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Fig. 4. — Datapoints show stacked flux densities (fx) of LAEs lack- 
ing individual detections in the SIMPLE IRAC images, with la error bars. 
The thick solid curve gives the best-fit model described in the text, which 
is a sum of a young coinponent (dashed curve) and an old component (thin 
solid curve). Squares show fluxes predicted by the best-flt model, which has 
xVd-o.f. = 14.6/12. 

i.e., a nearly constant star formation rate. Although we did 
not include our narrow -band photometry in the SED analysis, 
th e median LAE rest-fra me equivalent width of 60A found 
by IGronwall et al.l (l2007h is consistent with that expec ted for 
normal stellar populations in this age range (Finkelst ein et alj 
l2007h . The age of the old population is not well constrained 
but has a best fit of 2 Gyr (the age of the universe at z = 3. 1). 
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Fig. 5. — Constraints on age of the young stellar population versus its 
inass fraction. 

The mass fraction formed in the current starburst is not 
well constrained, and models where all of the stellar mass 
was produced in a current burst of star fo rmation of age 60 
to 350 Myr are allowed. Indeed, L ai et al.l d2007b) performed 
a one-component SED fit with .Marastonl (1200 5) models and 
found a best-fit age of 100 Myr, t=250 Myr, E(B-V)=0, and 
M» = 3 X 1O**M0. Our two-component best-fit is preferred to 
this, even accounting for the two extra degrees of freedom, but 
a single "r-model" population is not ruled out at 95% confi- 
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dence (see Fig. 5). 



5. DISCUSSION 



In CDM cosmology, galaxy formation is an ongoing pro- 
cess caused by merging of lower-mass dark matter halos, 
which may already possess stars. Finding stellar population 
ages of < 100 Myr is interesting. Our analysis of halo merger 
trees from the Milli-Millenium simulation ("Spring el et al.l 
l2005h found the median age of dark matter halos with M> 
IO'^-^Mq at z = 3.1 (defined as the age since half of the dark 
matter mass was accumulated) to be ^600 Myr, with only 
<10% of halos younger than 100 Myr. If repeated LAE 
phases occur, the mean halo occupation of ~ 1-10% can be 
interpreted as a "duty cycle" telling us what fraction of each 
halo's lifetime is spent in the early phases of starbursts be- 
fore significant dust is generated, and the population-averaged 
young age of ^^20 Myr would imply that this phase typically 
lasts ~40 Myr. Alternatively, if all dark matter halos expe- 
rience a single LAE phase shortly after their "formation" in 
a major merger, the mean halo occupation implies that LAEs 
will only be found in the youngest 1-10% of halos, which is 
barely consistent with their single-population best-fit age of 
100 Myr. If LAEs represent a subset of dark matter halos 
selected to have ages less than 100 Myr, their observed clus- 
tering may underestir nate their dark matter halo masses by up 
to a factor of two (see lGao &~Whit&.2007i) . 
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Fig. 6. — Tracks show the evolution of bias with redshift calculated using 
the no-merging model. The filled circle shows our result for the bias of LAEs 
at z = 3.1. Previous results at high-redshift are shown for L AEs at z = 4.5 
and z = 4.86 (stars, Kovac et al. 2007 and Ouchi et al. 2003, respectively), 
LBGs at z ~ 4 (stars, Ouchi et al. 2004), K-selected galaxies (diamonds. 



lOuadri et al.l2007l) . bright LBGs at z ~ 3 (triangle, Lee et al. 2006), and BM, 
BX, an d LBG galaxies (asterisks, Adelberger et al. 2005a, Adelberger et al. 
|2005b'). Local galaxy clustering is shown for SPSS galaxies (ope n circles, 
IZehavi et al. 2005) and for rich clusters fcross. iBahcall et ani2003ft . K-band 
limits are in Vega magnitudes. 

Figure |6] shows the reported bias values for LAEs to be 
lower than those of oth er z > 3 galaxy p opulations (bias 
values determined as in lOuadri et alj 120071) . The expected 



evolution of bias is show n for the "no-merging" model dFrvl 
Il996t IWhite et alJl2007h . A realistic amount of merging will 
cause the bias to drop somewhat faster, so the plotted trajec- 
tories provide an upper limit on the bias factor of a given 
point at lower redshifts. This shows that typical z = 3.1 
LAEs will evolve into galaxies of at most a few times L* 
at z = 0. The bias values imply that LAEs at z = 3.1 might 
evolve into the subset of BX galaxies at z — 2.2 dimmer 
than A" = 21.5, vvhich al so show relatively weak clustering 
(Adelberge r et alj|2005 bf). The K>21.5 BX galaxies have 
average M* = 1.5 x lO'^M©, so the z= 3.1 LAEs would need 
to form stars at an average rate of 14 MQyr"' over the inter- 
vening Gyr. This could be achieved with a constant specific 
SFR and no merging or with a constant SFR and ~2 major 
mergers. The only previous measur ements of LA E cluster- 
ing in u nbiased fields are a t z = 4.5 (iKovac e t al.' 2007) and 
z = 4.86 jOuchi et al.ll2003h . and this earlier LAE population 
appears to have significantly stronger clustering, consistent 
with possibly evolving into typical Lyman break galaxies at 

z~3. 

The models of iLe Delliou et all ( l2006l) predict stellar and 
dark matter masses and star formation rates for z = 3 LAEs 
within a factor of two of our results, despite assuming a 
top-heavy IMF and a very low escape fraction f^sc = 0.02 
that appears inconsistent with the observed lack of du st (see 
JK,objiyashi et al. 2007, for an alternative approach). Ma o et all 
JIOOTb used t he stel lar mass of 5 x 10 M© observed by 
iGawiser et al.l (l2006al) to predict LAE dark matter masses of 
10'° <M< IO^Mq, in the lower end of our allowed range. 
The stellar ages of ^ 20 Myr preferred by the two-population 
fit are noticeably lower than the maximu m values of 10 
to 500 My r predicted by these auth ors, iMori & Umemural 
(120061) . and lHaiman&SpaansI (Il999l) . but the ages of 60 to 
350 Myr preferred for the case of no evolved stars would be 
compatible. 

None of the current mod els a nd numeri cal simula- 
tions of LAEs (see also Thommes & Meisenheimer 2005|; 
iRazoumov & Sommer-Lars en 2006; Tasitsiomi,200 6) predict 
their present-day descendants. Nonetheless, the evolution of 
a significant fraction of z = 3.1 LAEs into z = L* galaxies 
with dark matter m ass Mom — 2 x lO'^M© and stellar mass 
M, ~ 4 X 10'°Mq (IIchikawaetalJ l2007) appears reasonable. 
Fig. |7] shows the histogram of present-day masses of dark 
matter halos in the Milli-Millenium merger trees that have 
progenitors with M> 5 x lO'^M© at z = 3.1. The median 
present-day halo mass is 1.2 x lO'^M©, and this would in- 
crease if LAEs found in sub-halos of massive z = 3.1 halos 
were included. Li et alJ (I2b07i) predict that the main progen- 
itor of a present-day L* galaxy had a dark matter mass of 
~ 1O"M0 at z = 3 and that these galaxies experienced sev- 
eral major mergers at 1.5 < z < 7. To form an L* galaxy at 
z = 0, several LAEs could merge while experiencing a mild 
reduction in average SFR, with accretion of lower-mass dark 
matter halos through minor mergers providing most of the fi- 
nal dark matter mass. However, the halo mass distribution of 
z = descendants in Fig|2]is very broad, with 25th and 75th 
percentile values of 2.9 x IO^'Mq and 7.6 x IO'^Mq. While 
z = L* galaxies like the Milky Way are roughly the median 
descendants of z = 3.1 LAEs, the descendant halos include a 
wide range from dwarf galaxies to rich galaxy groups. 

The typical LAE stellar mass at z = 3.1 is low er than that 
of any other studied high-redshift population (see Reddy et alJ 
,2006 ) but is close to that of dim (/ < 26 .3) Lyman break galax- 
ies (LBGs) at z ~ 5 (IVerma et al.ll2007h . LAEs at z = 3. 1 have 
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Fig. 7. — Histogram of dark matter halo masses of present-day descendants 
of halos with M> 5 X lO'^M© at z = 3.1. The dashed lines show the median 
halo mass of 1 .2 X IO'^Mq and the 25th and 75th percentile values of 2.9 X 
1O"M0 and 7.6 X lO'^M©. 

much lower star formation rate, stellar age, stellar mass, dark 
matter halo mass, and dust extinction than the ~ 3OM0yr"' , ~ 
500 Myr, - 2 x 10'"Mf7) ^ 3 x 10"Mo, Ay 1 LBG popula- 
tion at z - 3 (/? < 25.5. IShaplev et al.ll2001l:lAdelberger etal] 
|2005b ) or the - lOOMoyr"', ~ 2 Gyr, ~ lO'^Mp,, ^ lO'^Mfr, 
Ay c:i 2.5 Distant Re d Galaxy (DRG) p opulation (IWebb et alj 
120061 "Forster Schreiber et an i2004 : ijuadri et al.||2007l). The 
high-redshift Sub-Millimeter Galaxies {Chapman et al. 2003) 
appear to be the most massive and dusty, with the highest 
SFR. LAEs may represent the beginning of an evolutionary 
sequence where galaxies gradually become more massive and 
dusty due to mergers and star formation, but most LAEs at 
z = 3. 1 will never reach the DRG stage since DRG stellar and 
dark matter masses are akeady greater than those of present- 
day L* galaxies. 

The Damped Lya Absorption systems (DLAs, ' Wolfe et aP 
l2005h are another high-redshift population that probes the 
faint end of the luminosity function. The dark matter halo 
masses of DLAs at z 3 were determined by ICooke et al.l 



(!2_006!) to lie in the range lO'* <M< lO'^M© i.e., I3<b<4, 
which overlaps with the range of both L* and super-L* pro- 
genitors in Fig. |6] At least half of the DL As appear to have on- 
going star formation ( Wol fe et ani2004l) and two of the three 
DLAs detected in emission were seen in Lya. Further study 
is needed to determine the relationship between DLAs and 
LAEs. 

The observed properties of LAEs at z = 3.1 make them 
the most promising candidates to be high-redshift progenitors 
of present-day L* galaxies like the Milky Way. Our results 
suggest that LAEs are observed during the early phases of a 
burst of star formation, perhaps caused by a major merger of 
smaller dark matter halos. The input halos appear to have al- 
ready contained stars, accounting for the evolved stellar popu- 
lation that appears to contribute most of the LAE stellar mass, 
although starburst-only models are also allowed. It is clear 
that not all progenitors of L* galaxies were LAEs at z = 3.L 
The comoving number density of our sam ple of LA Es is a 
factor of 15 less than 0* for local galaxies (iLin et al.l ll996). 
plus we expect several high-redshift halos to merge into a sin- 
gle galaxy today. It remains possible that all progenitors of 
present-day galaxies experienced an LAE phase at some red- 
shift. Clustering and SED studies of LAEs at various redshifts 
are needed to assess the validity of this hypothesis. 
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